Low-Metallicity Star Formation: From the First Stars to Dwarf 
Galaxies 

Proceedings IAU Symposium No. 255, 2008 © 2008 International Astronomical Union 

L.K. Hunt, S. Madden & R. Schneider, eds. DOI: 00.0000/X000000000000000X 

Supernova Nucleosynthesis in the Early 

Universe 

Nozomu Tominaga 1 , Hideyuki Umeda 2 , Keiichi Maeda 3 , 
Ken'ichi Nomoto 3 ' 2 , and Nobuyuki Iwamoto 4 

1 Optical and Infrared Astronomy Division, National Astronomical Observatory, Mitaka, 

Tokyo, Japan 
email: nozomu . tominaga@nao .ac.jp 

2 Department of Astronomy, School of Science, University of Tokyo, Bunkyo, Tokyo, Japan 

3 Institute for the Physics and Mathematics of the Universe, University of Tokyo, Kashiwa, 

Chiba, Japan 

4 Nuclear Data Center, Nuclear Science and Engineering Directorate, Japan Atomic Energy 

Agency, Tokai, Ibaraki, Japan 

Abstract. The first metal enrichment in the universe was made by supernova (SN) explosions 
of population (Pop) III stars. The trace remains in abundance patterns of extremely metal- 
poor (EMP) stars. We investigate the properties of nucleosynthesis in Pop III SNe by means 
of comparing their yields with the abundance patterns of the EMP stars. We focus on (1) jet- 
induced SNe with various energy deposition rates [E dep = (0.3 — 1500) x 10 51 ergs s _1 ], and (2) 
SNe of stars with various main-sequence masses (M ms = 13 — 50Mq) and explosion energies 
[E — (1 — 40) x 10 51 ergs]. The varieties of Pop III SNe can explain varieties of the EMP stars: 
(1) higher [C/Fe] for lower [Fe/H] and (2) trends of abundance ratios [X/Fe] against [Fe/H]. 
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1. Introduction 

Long-duration 7-ray bursts (GRBs) have been found to be accompanied by luminous 
and energetic Type Ic supernovae [SNe Ic, called hypernovae (HNe)] (e.g. Galama et al. 
1998). Although the explosion mechanism is still under debate, photometric observations 
(a "jet break", e.g. Frail et al. 2001) and spectroscopic observations (a nebular spectrum, 
e.g. Maeda et al. 2002) indicate that they are aspherical explosions with jet(s). 

The aspherical explosions are also indirectly suggested from the abundance patterns 
of extremely metal-poor (EMP) stars with [Fe/H] < —30 The EMP stars are suggested 
to show nucleosynthesis yields of a single core-collapse SN (e.g. Beers & Christlieb 2005). 
Particularly, C-enhanced EMP (CEMP) stars have been well explained by the faint SNe 
with large fallback (Umeda & Nomoto 2005; Iwamoto et al. 2005; Nomoto et al. 2006; 
Tominaga et al. 2007b). On the other hand, some CEMP stars show enhancement of 
Co and Zn (e.g. Depagne et al. 2002) that requires explosive nucleosynthesis under high 
entropy. In a spherical model, however, a high entropy explosion is equivalent to a high 
energy explosion that inevitably synthesizes a large amount of 56 Ni, i.e., leads a bright 
SN (e.g. Woosley & Weaver 1995). This incompatibility will be solved if a faint SN is 
associated with a narrow jet within which a high entropy region is confined (Umeda & 
Nomoto 2005). 

f Here [A/B] = log 10 (iVA/iVB) — log 10 (iVA/iVB)©, where the subscript refers to the solar 
value and Na and Nb are the abundances of elements A and B, respectively. 
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Figure 1. (a) Schematic picture of the jet-induced explosion, (b) Density structure of the 40 
Mq Pop III star explosion model of i?de P ,5i = 15 at 1 sec after the start of the jet injection. 

2. Models 

We investigate the jet-induced explosions (e.g. Maeda & Nomoto 2003; Nagataki et al. 
2006) of 4OM Pop III stars (Umeda & Nomoto 2005; Tominaga et al. 2007b) using a 
two-dimensional special relativistic Eulerian hydrodynamic code (Tominaga 2009). 

The jets are injected at a radius i?o, corresponding to an enclosed mass of Mo, and the 
jet propagation is followed (Figs. [IJib). Since the explosion mechanism is unknown, the 
jets are treated parametrically with the following five parameters: energy deposition rate 
(£?dep), total deposited energy (E'dep)? initial half angle of the jets (0j e t)> initial Lorentz 
factor (rj et ), and the ratio of thermal to total deposited energies (/th)« 

In particular, we investigate the dependence of nucleosynthesis outcome on Edep fc> r 
a range of £?dep,5i = ^dep/10 51 ergs s _1 = 0.3 — 1500. The diversity of E^ep is consistent 
with the wide range of the observed isotropic equivalent 7-ray energies and timescales 
of GRBs (e.g. Amati et al. 2007). Variations of activities of the central engines, possibly 
corresponding to different rotational velocities or magnetic fields, may well produce the 
variation of i^dep- We expediently fix the other parameters as E'dep = 1.5 x 10 52 ergs, 
0jet = 15°, r jet = 100, / th = 10" 3 , and M = 1.4M (R ~ 900 km) in the models. 

The hydrodynamical calculations are followed until the homologously expanding struc- 
ture is reached (v oc r) and then the ejected mass elements are identified from their radial 
velocities. The nucleosynthesis calculations are performed as post-processing with ther- 
modynamic histories traced with marker particles that represent individual Lagrangian 
elements. In computing the jet composition, we assume that the jet initially has the 
composition of the accreted stellar materials. 



3. Jet-induced Supernovae 

3.1. Fallback 

Figure [2^ show "accreted" regions for models with i^dep,5i = 120 and 1.5 ergs s _1 , 
where the accreted mass elements initially located in the progenitor. The inner matter 
is ejected along the jet-axis but not along the equatorial plane. On the other hand, the 
outer matter is ejected even along the equatorial plane, since the lateral expansion of the 
shock terminates the infall as the shock reaches the equatorial plane. 

The remnant mass (M rem ) is larger for lower i^dep- This stems from the balance between 
the ram pressures of the injecting jet (Pj e t) and the infalling matter (Pf a ii). In order to 
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Figure 2. (a) Initial locations of the mass elements which are finally accreted for models with 
^de P ,5i = 120 (black) and with Ed ep ,5i = 1.5 (cyan), (b) Ejected Fe mass (solid line: explosive 
nucleosynthesis products, dotted line: the jet contribution) as a function of the energy deposition 
rate, (c) Dependence of abundance ratio [C/Fe] on the energy deposition rate. 

inject the jet, Pj et should overcome Pf a ii- Pjet is determined by Ro, i^dep, #jet> Tj et , and 
/th, thus being constant in time in the present models. On the other hand, Pf a n decreases 
with time, since the density of the outer materials decreases following the gravitational 
collapse (e.g. Fryer & Meszaros 2003). For lower Edep, ^jet is lower, so that the jet 
injection (Pj et > Pf a ii) is realized at a later time when the central remnant becomes more 
massive due to more infall. As a result, the accreted region and M rem are larger for lower 

^dep- 

A model with lower i^dep has larger M rem , higher [C/Fe], and smaller amount of Fe 
[M(Fe)] because of the larger amount of fallback (Figs.[2j3c, Tominaga et al. 2007a). The 
larger amount of fallback decreases the mass of the inner core relative to the mass of 
the outer layer. The fallback of the O layer also reduces M(Fe) because Fe is mainly 
synthesized explosively in the Si and O layers. The variation of E^ep in the jet-induced 
explosions predicts that the variation of [C/Fe] corresponds to that of M(Fe). 

3.2. Comparison with the spherical supernova model 

The calculations of the jet-induced explosions show that the ejection of the inner matter 
is compatible with the fallback of the outer matter (Fig. [2^). This is consistent with the 
two-dimensional illustration of the mixing-fallback model (Fig. [3fci) proposed by Umeda 
& Nomoto (2002). The mixing-fallback model has three parameters; initial mass cut 
[M cut (ini)], outer boundary of the mixing region [M mix (out)], and a fraction of matter 
ejected from the mixing region (/). The remnant mass is written as 

M rem = M cut (ini) + (1 - /)[M mix (out) - M cut (ini)]. (3.1) 

The three parameters relate to the hydrodynamical properties of the jet-induced explo- 
sion. 

A model with E dePj51 = 120 and M = 2.3M (R ~ 3 x 10 3 km) is compared with 
the spherical SN explosion adopting the mixing-fallback model. The inner boundary 
and the central remnant mass of the jet-induced explosion model are Mq = 2.3M and 
M rem = 8.1M . Its abundance pattern is well reproduced by the spherical SN model 
with the same main-sequence mass M ms = 4OM and the explosion energy E = 3 x 10 52 
ergs (Fig. Eb). The parameters of the mixing-fallback model are M cut (ini) = 2.3M , 
M mix (out) = 1O.8M and / = 0.19. The resultant M rem (= 9.2M ) is slightly larger 
than M rem of the jet-induce SN model. 
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Figure 3. (a) Two-dimensional illustration of the mixing- fallback model, (b) Comparison of the 
abundance patterns of the jet-induced SN model with i?de P ,5i = 120 and Mo = 2.3M© (solid 
line) and the mixing- fallback model (dashed line). 



There, however, are some elements showing differences. The differences stem from 
the high-entropy explosion due to the concentration of the energy injection (e.g. Maeda 
& Nomoto 2003). The enhancements of [Sc/Fe] and [Ti/Fe] improve agreements with 
the observations. Such thermodynamical feature of the jet-induced explosion model can- 
not be reproduced by the mixing-fallback model, while a "low-density" modification 
might mimic the high-entropy environment (e.g. Umeda & Nomoto 2005; Tominaga et al. 
2007b). 



4. Trends with Metallicity 

The abundance patterns of the EMP stars show certain trends of abundance ratios 
[X/Fe] with respect to [Fe/H] (Cayrel et al. 2004). To clarify the origin of the trends, 
SN models with M ms = 13 — 5OM are calculated. The explosion energies are set to be 
consistent with the observations of present SNe (e.g. Tanaka et al. 2008). 

The abundance ratios against [Fe/H] are compared with yields of individual SN models 
and the Salpeter's IMF-integrated yield (Fig. [4]). [Fe/H] of a next-generation star is 
determined by M(Fe) and a swept-up H mass. Since the swept-up H mass is almost 
proportional to E of the SN, [Fe/H] of the next-generation stars are determined by a 

. 6 / 7 l 



equation [Fe/H] = log 10 



- C (Thornton et al. 1998), where C is 
of the IMF- integrated abundance ratios 



M(Fe) // E V 
M© / \^10 51 ergs J 

assumed to be a constant value of 1.4. [Fe/H 
are assumed to be the same as SN models with E ~ 10 51 ergs ([Fe/H] ~ —2.6). 

HNe explode with E > 10 52 ergs and eject large amount of 56 Ni (> O.1M ), while 
normal SNe explode with E ~ 10 51 ergs and eject ~ 0.07Mq of 56 Ni. Therefore, according 
to the above equation, [Fe/H] of a next-generation star originated from a HN is lower than 
that of a next-generation star originated from a normal SN. The higher-energy explosion 
raises explosive nucleosynthesis under higher entropy and thus leads higher [Zn/Fe]. This 
accompaniment explains the observed trends of [Zn/Fe] against [Fe/H]. The trends of 
other elements are also reproduced by the variations of M ms and E (Fig. [4]). 



5. Conclusion 

We focus on two interesting properties observed in the abundance patterns of the metal- 
poor stars: (1) the higher [C/Fe] for lower [Fe/H] and (2) the trends of [X/Fe] against 
[Fe/H]. The variations of the metal-poor stars are explained by the variations of SNe 
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Figure 4. Comparison between the [X/Fe] trends of observed stars (e.g. Cayrel et al. 2004; 
Honda et al. 2004: cross) and SN models with the mixing-fallback model and applied the Y e and 
"low-density" modifications (Tominaga et al. 2007b: circles). 

that contribute the metal enrichment of the early universe. Especially, (1) the variation 
of the energy deposition rates explains the tendency of [C/Fe] against [Fe/H] and (2) the 
variations of M ms and E explain the trends of [X/Fe] against [Fe/H]. We propose that 
the abundance patterns of the metal-poor stars will provide additional constraints on 
the explosion mechanism of GRBs and SNe other than the direct observations of present 
GRBs and SNe. 
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